Heavy-quark production in proton-nucleus collisions at the LHC 
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A sizable rate of events, with several pairs of heavy-quarks produced contemporarily by 
multiple parton interactions, may be expected at very high energies as a consequence of the 
large parton luminosities. The production rates are enhanced in hadron-nucleus reactions, 
which may represent a convenient tool to study the phenomenon. We compare the different 
contributions to bbbb, cccc and bbcc production due to single and double parton scatterings, 
in collisions of protons with nuclei at the CERN-LHC. 
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(N '. 1. INTRODUCTION 

£> ' The large rates of production of heavy quarks, expected at high energies, may lead to a sizable 
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number of events, at the CERN-LHC, containing different pairs of heavy quarks, generated 
contemporarily by independent partonic collisions. An inclusive cross section of the order of 10 
fib may in fact be foreseen for a double parton collision process, with two bb pairs produced 



■"sj - . in a interactions at 14TeV[l], while the cross section to produce two cc pairs may be of 



the order of one mb, the contribution of single parton collisions to the processes being one 
order of magnitude smaller. All production rates are significantly enhanced in proton-nucleus 
collisions, which may offer considerable advantages for studying multiparton collision[2]. Given 
^| the large rates expected, the production of multiple pairs of heavy-quarks should hence represent 

a convenient process to study multiparton interactions in pA collisions at the LHC. On the other 
hand the mechanism of heavy quarks production is not yet understood satisfactorily also in the 
^ ■ simplest case of nucleon-nucleon collisions, the effects of higher order corrections in a s being still 

a matter of debate. A comprehensive description, of the much more structured process of heavy 
quarks production in hadron-nucleus interactions, might hence be approached after gaining a 
deeper understanding of the short scale parton-level dynamics of the process. On the other 
hand a significant feature of higher order corrections in a s is that, for a limited set of physical 
observables, the whole effect of higher orders reduces to an approximate rescaling of a lowest 
order calculation of heavy quarks production in perturbative QCD [3]. Some of the features of 
the process are therefore effectively described by the simplest parton level dynamics at the lowest 
order in a Sl which lets one speculate that a similar property might hold also for a much more 
complex process, as multiple production of heavy quarks. 
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Taking this optimistic point of view we will attempt, in the present note, to obtain indications 
on some properties of bbbb, cccc and bbcc production in pA interactions at the LHC, by considering 
the contributions of the two different interaction mechanisms, the connected 2^4 and the 
disconnected (2 — > 2) 2 parton processes, which will be effectively described by the lowest order 
diagrams in perturbative QCD, while keeping higher order corrections into account by a simple 
overall rescaling. Obviously, following this philosophy, all considerations are necessarily limited 
to the restricted class of physical observables falling in the category above. 

2. PRODUCTION IN PROTON-NUCLEUS COLLISIONS 

Quite in general [4, 5] with the only assumption of factorization of the hard component of the 
interaction, the expression of the double parton scattering cross section to produce two pairs of 
heavy-quarks is given by 

a (p,A) = ~~2 / ^p( Xi ' 1 x i' s ij)®( x ii x i)&( x j > x> j)T(N,A)( x 'ji x> j> Sij)dxidXidxjdXjd 2 Sij, (1) 

ij 

where the index N refers to nucleon and A to a nucleus, while the indices i,j to the different 
kinds of partons that annihilate to produce a given qq pair and the factor m/2 is a consequence 
of the symmetry of the expression for exchanging i and j; specifically m = 1 for indistinguishable 
parton processes and m = 2 for distinguishable parton processes. The interaction region of a 
hard process is very small as compared to the hadron scale, hence in the case of a double parton 
collision the two elementary interactions are well localized in transverse space, within the two 
overlapping hadrons. The rate of events where two hard collisions take place contemporarily, in a 
given inelastic hadron-hadron event, depends therefore on the typical transverse distance between 
the partons of the pairs undergoing the multiple processes. A main reason of interest is hence that 
double parton scatterings may provide information on the typical transverse separation between 
pairs of partons in the hadron structure. Indeed the non-perturbative input of a double parton 
collision is the two-body parton distribution function T(xi, X2, ^1,2), which depends not only on 
the fractional momenta x± t 2, but also on the relative distance in transverse space si ; 2, besides 
(although not written explicitly to simplify the notation) on the scales of the two interactions and 
on the different kinds of partons involved. As a consequence the double parton scattering cross 
section is characterized by a linear dependence on dimensional scale factors, which are related 
directly with the typical transverse distances between the various pairs of partons, contributing 
to the double scattering process under consideration. 

The cross section is simplest when the target is a nucleon and partons are not correlated in 
fractional momenta, which may be a sensible approximation in the limit of small x. In such a 
case the two-body parton distribution may be factorized as T p (xi, Xj] Sij) = G(xi)G(xj)F(sij), 
where G(x) is the usual one-body parton distribution and F(s) a function normalized to 1 and 
representing the parton pair density in transverse space. With this simplifying assumptions the 
inclusive cross section to produce two pairs of heavy quarks is written as [6] 

ij 
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where Oi Oj represent the inclusive cross sections to produce a qq pair in a hadronic collision, with 
the indexes i, j labelling a definite parton process. The factors Q lJ have dimension an inverse 
cross section and result from integrating the products of the two-body parton distributions in 
transverse space. In this simplified case the factors Q tJ provide a direct measure of the different 
average transverse distances between different pairs of partons in the hadron structure [6, 7]. 

The cross section has a more elaborate structure in the case of a nuclear target. The most 
suitable conditions to study the phenomenon are those where the nuclear distributions are addi- 
tive in the nucleon parton distributions. In such a case one may express the nuclear parton pair 
density, Ta(x'j, x'p Sij), as the sum of two well defined contributions, where the two partons are 
originated by either one or by two different parent nucleons: 



^ A(%ii Xj', Sij) — r^(xj, Xj ; Sij) ^ + r^(xj, Xj ; s 



(3) 



and correspondingly = cr^ 1 1 + cr^, 1 2 • The two terms r^ji^ are related to the nuclear nucleon's 
density. Introducing the transverse parton coordinates B±?rf, where B is the impact parameter 
of the hadron-nucleus collision, one may write 



r a (%i 1 Xj ; s t j t 



1,2 



(4) 



where 7a|i,2 are given by 



lA[x'i,x'j-B 



D J2_ 
2 ' 



B 



lA (x' t ,x' f ,B+ S f,B-f 



r N (x'i,x'j;sij)T(B) 



G N {4)G N {x'j)T(B+ S -f)T(B- S -f) 



(5) 



with T(B) is the nuclear thickness function, normalized to the atomic mass number A and Gn 
nuclear parton distributions divided by the atomic mass number. 

In the simplest additive case, the first term in Eq.(3) obviously gives only a rescaling of the 
double parton distribution of a isolated nucleon 



TA^x'j-Sij^T^x'j-Sij) j d 2 BT(B) 



(6) 



and the resulting contribution to the double parton scattering cross section is the same as in a 
nucleon-nucleon interaction, apart from an enhancement factor due to the nuclear flux, which is 
given by the value of the atomic mass number A: 



° D A 



Aa D 



N- 



(7) 



The c;^ 1 2 term has more structure. In this case the integration on the relative transverse distance 
between the partons of the interacting pairs, s^, involves both the projectile and two different 
target nucleons: 



j ds ij T p {x i ,Xj-s ij )T(B+ S -ipjT(B 



(8) 
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As one may notice the expression (8) depends on two very different scales, the hadron radius r p 
and the nuclear radius Ra- A usual approximation in pA interactions is to consider the limit 
r p <C Ra, where one may use the approximation 

r(B±|)^T(B), (9) 

which allows one to decouple the integrations on and on B. One hence obtains: 



where 



- y~] / Gp(xi,Xj)a(xi,x' i }a(xj,x'j)GN(x' i )GN(x'j)dxidx' i dxjdx'j I d 2 BT 2 (B), (10) 



Gp(xi,Xj*) — / d SijJ^p^Xi, Xj] Sjj) . 



(11) 



Remarkably the two terms a® 



and cr^ 



have very different properties. In fact, while the 



correct dimensionality of 



is provided by transverse scale factors related to the nucleon scale, 



eqs.(2, 7), the analogous dimensional factor is provided in a® 



by the nuclear thickness function, 



which is at the second power in Eq.(10), being two the target nucleons involved in the interaction. 
As pointed out in ref. [2], while on general grounds &f p A) depends both on the longitudinal 

and transverse parton correlations, the a® term depends solely on the longitudinal momentum 

fractions that, when the a®] term is isolated one has the capability of measuring the 

longitudinal and, a fortiori, also the transverse parton correlations of the hadron structure in a 
model independent way 

Although the two contributions may be defined in a more general case, the separation of the 



cross section in the two terms a® 



and a® 



is most useful in the regime of additivity of the 



nuclear structure functions, which may not be a bad approximation for a sizable part of the 
kinematical regime of heavy-quarks production at the LHC. In the case of a central calorimeter 
with the acceptance of the ALICE detector (\t]\ < 0.9), the average value of momentum fraction 
of the initial state partons, in a pp — > bbbb process, is (x) ~ 6 x 10 -3 while for pp — ► cccc one 
finds (x) « 2 x 10~ 3 . With a cut of 5 GeV in the transverse momenta of the b, c quarks one 
obtains (x) xs 10~ 2 and (x) « 8 x 10~ 3 respectively, while a cut of 20 GeV gives (x) « 3 x 10 -2 . 
In the case of 2.5 < rj < 4 and without any cut in p t , the average values of momentum fraction 
are (x) « 3 x 10 -2 and (x) rs 9 x 10~ 3 [8]. Deviations from additivity at low x are less than 
10% for x > 2 x 10~ 2 [10] and, although increasing with the atomic mass number, non additive 
corrections should hence be at most a 20% effect, in above the kinematical regimes. 

As mentioned in the introduction, heavy quarks production is characterized by a non trivial 
dynamics, in such a way that also next-to-leading corrections to the lowest order term in pertur- 
bative QCD are not sufficient for an exhaustive description of the inclusive spectra, which most 
likely need an infinite resummation to be evaluated. After comparing the results of different 
approaches to heavy quarks production (as NLO QCD and fci-factorization) one nevertheless 
finds that, in a few cases, the whole effect of higher orders is to a large extent just an approxi- 
mate rescaling of the results obtained by a lowest order evaluation in perturbative QCD. When 
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limiting the discussion to an accordingly restricted set of physical observables, the whole effect 
of higher order corrections to heavy quark production is hence summarized by a single number, 
the value of the so-called K-factor defined as: 

K=«*iL (12) 

where cr(qq) is the inclusive cross section for qq production and (Tlo(qq) the result of the low- 
est order calculation in pQCD. By evaluating, with the fcj-factorization approach, rapidity and 
pseudorapidity distributions of bb and of cc production, in pp collisions at the center-of-mass 
energy of y/s = 5.5 TeV, within \rj\ < 0.9 and 2.5 < rj < 4, with different choices of heavy-quark 
masses and of the factorization and renormalization scales, one finds a result not incompatible 
with a lowest order calculation in perturbative QCD rescaled by the factors K shown in Tables 
[I,II,III]. 

To obtain the values of the K-factor in Table [I], in the evaluation of the cross section with 
the /^-factorization approach, we have set the factorization scale equal to the invariant mass of 
the parton process s and the renormalization scale equal to the virtuality of the initial state 
gluons, while, for the cross section at the lowest order in pQCD, we used as a scale factor the 
transverse mass of the produced quarks. In Tables [II, III] we rather used in both cross sections 
the average of the squared transverse masses of the produced quarks and the heavy-quark mass 
as scale factors [8]. 

With a fixed choice of factorization and renormalization scales, one obtains a substantial 
increase of the value of the K factor when decreasing the value of the heavy quark mass, the 
increase being larger for cc than for bb production. With our different choices for the scale factors 
we obtain variations of the K factor within 4% for bb and 15% for cc production. 

To discuss qqqq production in pA collisions, while remaining in a kinematical regime where 
non additive corrections to the nuclear structure functions are not a major effect, we have limited 
all considerations to physical observables, where higher orders may be taken into account by a 
simple rescaling of the lowest order calculation. For qqqq production, where only results of three 
level calculations are up to now available, we have further assumed that the value of the K- 
factor in the 2^4 parton process is the same as in the 2^2 process [1, 9]. We have hence 
evaluated the various contributions to the cross section in the case of a central, \r]\ < 0.9, and 
of a forward calorimeter, 2.5 < rj < 4, where the approximate expression of the cross section 
in Eq.(2) may not be an unreasonable approximation. As the process is dominated by gluon 
fusion, the expression of the cross section in Eq.(2) may be limited to a single term only. For the 
corresponding dimensional scale factor we used the value reported by the CDF measurement of 
double parton collisions [11, 12] and to evaluate the cross section at the lowest order in pQCD 
we used the MRS99 parton distributions [13], with factorization and renormalization scale equal 
to the transverse mass of the produced quarks. The cross sections of the 2^4 processes have 
been evaluated generating the matrix elements of the partonic amplitudes with MadGraph [14] 
and HELAS [15]. For the mass of the bottom quark we used the central values = 4.6 and 
m c = 1.4 GeV. The multi-dimensional integrations have been performed by VEGAS [16] and the 
lowest order pQCD cross sections have been finally multiplied by the K factors of Table [I]. 
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3. RESULTS 

A major result of the present analysis is that the effects induced by the presence of the 
nucleonic degrees of freedom, in double parton scatterings with a nuclear target, cannot be 
reduced to the simple shadowing corrections of the nuclear parton structure functions, whose 
effect is to decrease the cross section as a function of A. In the case of double parton collisions, 
the main effect of the nuclear structure is represented by the presence of the term in the 
cross section, which scales with a different power of A as compared to the single scattering 
contribution, producing an additive correction to the cross section. 

To emphasize the resulting "anomalous" dependence of the double parton scattering cross sec- 
tion, as a function of the atomic mass number, we have plotted in fig.[l] the ratio at? / (erf + a s ), 
as a function of A. The ratio represents the contribution to the cross section of the processes 
where two different nuclear target nucleons are involved in the interaction, scaled with the con- 
tribution where only a single target nucleon is involved. The dependence on the atomic mass 
number of the latter terms of the cross section, namely the single (a s ) and the double (erf) 
parton scattering terms against a single nucleon in the nucleus, is the same of all hard processes 
usually considered, where nuclear effects may be wholly absorbed in the shadowing corrections 
to the nuclear structure functions. The contribution to the cross section of the erf term is, on 
the contrary, "anomalous", involving two different target nucleons in the interaction process. The 
ratio above hence represents the relative weights of the "anomalous" to the "usual" contributions 
to the double parton scattering cross section on a nuclear target. 

The plots in Fig. [2] show the rapidity distribution in a forward calorimeter (2.5 < rj < 4) of a 
heavy quark produced in an event with bbbb, cccc and bbec respectively: one-nucleon (dashed 
lines) and two- nucleon contributions (continuous lines) in the case of a heavy, Au, (higher curves) 
and of a light nucleus, O, (lower curves). In each case considered continuous and dashed curves 
differ essentially only by a rescaling, showing that the effect of the single parton scattering term 
(the 2^4 parton process) is negligible in this kinematical regime. The different contributions to 
the cross section for bbbb, cccc, bbec production, due to interactions with a single (dashed line) or 
with two different target nucleons (dot line), are shown in Fig. [3,4,5] as a function of the atomic 
mass number in the case of a central and of a forward calorimeter, with two different choices 
of cuts (pt = 0, lOGeV) on the outgoing quark-transverse momenta. The continuous line is the 
sum of the two contributions. As one may see, by introducing a cutoff in pt one enhances the 
single scattering contribution and the relative role of af is decreased. 

To have some indication on the overall uncertainty of our estimates we have plotted in Fig. [6] 
the range of values obtained for the integrated cross sections, for bbbb, cccc and bbec production, 
within \rj\ < .9, as a function of the atomic mass number, when making the different choices 
described above. In each panel we report the different choices corresponding to the extreme 
values of the cross section. The indication obtained in this way is that the overall cross section 
is roughly determined within a factor three. 

Summarizing the large size of the cross section of heavy-quark production in hadron-nucleus 
collisions at the LHC (one may expect values of the order of 5 — 10 mb for charm production) 
suggests that the production of multiple pairs of heavy quarks is fairly typical at high energies, 
hence representing a convenient channel to study multiple parton interactions. A rather direct 
feature, which is a simplest prediction and then a test of the interaction mechanism described 



7 



above, is the "anomalous" dependence on A. The effects induced by the presence of the nucleonic 
degrees of freedom in the nuclear structure are in fact not limited, in this case, to the shad- 
owing corrections to the nuclear structure functions usually considered, which cause a limited 
decrease (not larger than 20%, in the kinematical regime considered here) of the cross section for 
a hard interaction in hadron-nucleus collisions. When considering double parton scatterings, all 
nuclear effects can be exhausted in the shadowing corrections to the nuclear structure functions 
in the <J^\i term only. The dominant effect of the nuclear structure is on the contrary due to 
the presence of the cr^ (2 term in the cross section, which scales with a different power of A as 
compared to single scattering term, giving rise to a sizably larger correction, with opposite sign 
as compared to the shadowing correction, namely to an increase of the cross section, which may 
become larger than 100% for a heavy nucleus. 
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fi F — s /x^ — q 2 (gluon virtuality) 


m b (GeV) 


K b = <T(bb)/cr L o(bb) 


m c (GeV) 


K c = o-(cc)/o- LO (cc) 


4.25 


5.9 


1.2 


7.0 


4.5 


5.7 


1.4 


6.6 


4.75 


5.5 


1.6 


6.1 



Table I: JiT factor values for different choices of the mass of the heavy quarks. In the k t factorization 
approach the cross section has been evaluated by using s as factorization scale and the gluon virtuality 
as renormalization scale; when evaluating the cross section at the lowest order in pQCD the scales have 
been set equal to the transverse mass of the heavy quarks [8] 



Hp = »R = Mo = ( m l Q + m lo)/ 2 


m b (GeV) 


K b = a(bb)/cr L o(bb) 


ra c (GeV) 


K c — cr(cc)/cr LO (cc) 


4.25 


6.1 


1.2 


6.2 


4.5 


5.9 


1.4 


5.9 


4.75 


5.8 


1.6 


5.6 



Table II: K factor values for different choices of the masses of the heavy quarks. In the calculation of the 
cross sections, factorization and renormalization scales have been set equal to the average of the squared 
transverse masses of the heavy quarks[8] 



/4 = = m h 


m b (GeV) 


K b = a(bb)/cT L o(bb) 


m c (GeV) 


K c — a-(cc)/cr L o(cc) 


4.25 


5.7 


1.2 


8.1 


4.5 


5.6 


1.4 


7.9 


4.75 


5.4 


1.6 


7.3 



Table III: K factor values for different choices of the masses of the heavy quarks. In the calculation of 
the cross sections, factorization and renormalization scales have been set equal to the mass of the heavy 
quark. [8] 



pA -> bbbb Vs = 5.5 GeV pA cccc Vs = 5.5 GeV 



pA -> bbcc Vs = 5.5 GeV 



- : 1 11 1 < 0.9 p cut = GeV 

■ : It] I < 0.9 p cul = 10 

. :2.5 <T)< 4 p cut = GeV 

. :2.5 <11< 4 p cut = 10 GeV 

: 4.5 GeV 




j I i I i I i I i I i i o 

40 80 120 160 200 240 
A 



- : h 1 < 0.9 p ci 


, = GeV 


. 1 T| | < v.y p ci 


_ 1Q Q e V — 




= GeV 


~ :2.5<ri<4p 


= 10 GeV " 


~ m = 1.4 GeV 

c 




«< 

_ AT 








T 


1 i 1 i 



H | < 0.9 p ci 
11 1 < 0.9 p 



= 0GeV 
= 10 GeV 



. :2.5 <T1< 4 p cut = GeV 
. :2.5<Ti<4p ci]t =10GeV 

= 4.5 GeV 
1.4 GeV 




40 80 120 160 200 240 
A 





40 80 120 160 200 240 
A 



Figure 1: Relative weights of the terms with "anomalous" and "usual" ^4-dependence in the double scat- 
tering cross section for bbbb, cccc, bbcc production. 
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Figure 2: Pseudorapidity distribution in a forward calorimeter of a heavy quark produced in events with 
bbbb, cccc and bbcc: one-nucleon (dashed lines) and two-nucleon contributions (continuous lines) in the 
case of a heavy (higher curves) and of a light nucleus (lower curves) . 
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Figure 3: Different contributions to the cross section of bbbb production in a central and in a forward 
calorimeter as a function of A. Cross sections without any cut in p t (left figures) and after applying a cut 
of 10 GeV in the transverse momenta of each produced heavy-quark (right figure): one-nucleon (dashed 
lines), tv 




Figure 4: Different contributions to the cross section of cccc production in a central and in a forward 
calorimeter as a function of A. Cross sections without any cut in p t (left figures) and after applying a cut 
of 10 GeV in the transverse momenta of each produced heavy-quark (right figure): one-nucleon (dashed 
lines), two nucleons contribution (dotted lines) and total cross section (continuous lines). 
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Figure 5: Different contributions to the cross section of bbcc production in a central and in a forward 
calorimeter as a function of A. Cross sections without any cut in p t (left figures) and after applying a cut 
of 10 GeV in the transverse momenta of each produced heavy-quark (right figure): one-nucleon (dashed 
lines), two nucleons contribution (dotted lines) and total cross section (continuous lines). 
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Figure 6: Production cross section of bbbb, cccc and bbcc in a central calorimeter for different choices of 
quark-masses and factorization scales. Lower and upper curves, correspond respectively to the lower and 
upper values of the heavy quark masses considered in our calculation. 



